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Abstract. We study the motion of a particle moving on a two-dimensional 
honeycomb lattice, whose sites are randomly occupied by either right or left 
rotators, which rotate the particle’s velocity to its right or left, according to 
deterministic rules. In the model we consider, the scatterers are each initially 
oriented to the right with probability p E [0,1]. This is done independently, so 
that the initial configuration of scatterers, which forms the medium through 
which the particle moves, are both independent and identically distributed. 

For p E (0,1), we show that as the particle moves through the lattice, it cre¬ 
ates a number of reflecting structures. These structures ultimately limit the 
particle’s motion, causing it to have a periodic trajectory. As p approaches 
either 0 or 1, and the medium becomes increasingly homogenous, the particle’s 
dynamics undergoes a discontinuous transition from this self-limiting, periodic 
motion to a self-avoiding motion, where the particle’s trajectory, away from its 
initial position, is a self-avoiding walk. Additionally, we show that the periodic 
dynamics observed for p E (0,1) can persist, even if the initial configuration of 
scatterers are not identically distributed. Furthermore, we show that if these 
orientations are not chosen independently, this can drastically change the par¬ 
ticle’s motion causing it to have a behavior that is nonperiodic. 

PACS numbers: 05.50-|-q, 02.10.Ox 


1. Introduction 


The particular system we consider in this paper is an example of a Lorentz lattice 
gas (LLG). In a LLG, a single particle moves along the bonds of a lattice, from 
lattice site to lattice site. When the particle arrives at a lattice site, it encounters 
a scatterer that modifies its motion according to a given scattering rule. 

The reason to study such systems is to understand the basic principles that un- 
derly dynamic processes such as diffusion, propagation, etc. [2-10 . For simplicity, 
the study of a particle’s motion on a lattice is a natural choice, since a lattice has 
both a discrete structure and a high degree of regularity. In such systems the type 
of scattering rules that have been investigated are physically motivated rules such 
as rotators, mirrors, 


etc. 


11-14 


In the case that there is a scatterer at each lattice site and each scatterer is fixed, 
i.e. is not affected by the particle’s motion, the problem of determining the particle’s 
motion through the lattice is related to problems in percolation theory [5,12 15 . 
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When the scatterers are not fixed, as is the case in this paper, and are affected 
by the particle, the particle’s motion is a much more dynamic process and has 


connections to problems in kinetic theory 16,17 . 

In this paper, we continue our investigation, begun in |^, of the particle’s motion 
on the regular two-dimensional honeycomb lattice, in which the lattice is fully 
occupied by flipping scatterers. The particular type of scatterers we consider here 
are flipping rotators, which rotate the particle’s velocity either to its left or its right 
by an angle of 0 = ±7r/3, depending on whether the scatterer is oriented to the left 
or the right, i.e. is a left or right scatterer, respectively. Furthermore, the scatterers 
flip or change orientation after scattering the particle, flipping either from right to 
left or from left to right, depending on their original orientation, respectively. 

This flipping motion of the scatterers has a number of consequences. Viewing 
the lattice and its scatterers as a medium through which the particle moves, the 
fact that the particle can change (flip) the scatterer’s orientation implies that there 
is an interplay or interaction between the particle and the medium. That is, the 
particle affects the medium, which then in turn, has an effect on the particle’s 
motion. It is worth mentioning that this type of system, which we consider in this 
paper, is likely one of the simplest systems in which there is an interaction between 
the particle and medium. 

Having fixed a lattice and a scattering rule, it remains to choose an initial con¬ 
figuration of the scatterers, in order to study the properties of the motion of the 
particle over the lattice. The initial configuration we consider, throughout the ma¬ 
jority of this paper, is the configuration in which each scatterer is a right rotator 
with probability p € [0,1]. The orientation of each scatterer is chosen independent 
from the others, so that the collection of their initial orientations forms a set of in¬ 
dependent identically distributed (i.i.d.) random variables. From a physical point 
of view, we interpret this as an assumption that the particle moves in a random 
medium. 

Although the initial configuration is random, the particle’s motion through the 
lattice is governed by deterministic scattering rules. In this setting, the particle’s 
motion is referred to as a deterministic walk in a random medium or environment 
That is, the first time the particle visits a lattice site, it is scattered to its 


18 


right by a right scatterer with probability p S [0, !]• However, on each subsequent 
visit to this same site, the particle is scattered (deterministically) in the opposite 
direction, from which it was previously scattered, since the scatterer has flipped its 
orientation after each visit. 

We find that, as long as p S (0,1), the particle’s motion will be periodic (see 
theorem 2.3). Because the random initial configuration in this model initially causes 
the particle to move randomly, this means that, the deterministic nature of the 
particle’s motion eventually overcomes this randomness. The reason this takes 
place is that, as the particle moves through the lattice, it creates a number of 
structures, which we refer to as reflecting structures. The initial effect of a reflecting 
structure on the particle is to cause the particle to reverse its trajectory back to its 
initial position. However, a second, and arguably more important effect of these 
structures, is to block the particle’s progress through the lattice. We find that a 
number of reflecting structures, acting together, can permanently limit the particle’s 
motion, causing its trajectory to become periodic (see theorem 3.7). Thus, we refer 
to this particular type of periodic motion as self-limiting. 
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Reflecting structures and their effect on the particie’s motion have been previ- 
ousiy observed in both the square 11 and the honeycomb lattice [^. In this paper, 
we can go further and describe the general topology of these structures. We also 
introduce two other structures and describe their topology and effect on the par¬ 
ticle. These new structures are, respectively, called semi-reflecting structures and 
reflector transforms. 

We show that as the particle moves through the lattice, it can either create, 
transform, or annihilate a reflecting structure. In order to become trapped or limited 
by a number of reflecting structures, we show that the particle must go through 
the process of both creating as well as transforming these structures. However, 
to avoid becoming trapped, the particle can also annihilate those reflectors it has 
either created or transformed. In this sense, the particle acts like an architect that 
can build, destroy, and remodel, i.e. transform, these structures. 

This competition between reflector formation, transformation, and annihilation, 
occurs for all p-values in the set (0,1). However, for p = 0,1 the situation is very 
different. Building on recent results [^, we show that for these two p-values, the 
particle can never create a reflecting structure (see corollary]^. Additionally, we 
observe that the lack of reflecting structures, causes the particle to have a non¬ 
periodic and therefore unbounded trajectory [^. Hence, as p approaches 0 or 1, i.e. 
as the random medium becomes increasingly homogenous, the particle undergoes 
a discontinuous dynamical transition from having a periodic, self-limiting motior0 
to having a very different non-periodic behavior. 

In the periodic, or self-limiting case, the particle becomes trapped between the 
reflecting structures it creates. Once the particle creates a reflecting structure, it 
will reverses its entire trajectory back to its initial position, going only through 
those lattice sites it has already visited. In contrast, for p = 0,1, the particle not 
only appears numerically to have an unbounded trajectory but, between returns to 
its initial position, the particle never visits any lattice site twice (see theorem 5.1). 
This self-avoiding motion is then vastly different from the periodic self-limiting 
motion observed for p = (0,1). 

In what seems paradoxical, in a random medium, i.e. for p S (0,1), the particle 
has a motion that is ordered in the sense that its trajectory is periodic and therefore 
regularly repeats itself. On the other hand, when p = 0,1 and the medium is 
homogeneous, and therefore much more ordered, the particle’s motion is much more 
irregular exhibiting features of deterministic chaos and of self-organized criticality 

01 ^ ■ 


To more fully understand how and why reflecting as well as semi-reflecting struc¬ 
tures appear in different random media, i.e. for different p € (0,1), we will also 
consider a generalization of our flipping rotator model, in which the scatterer’s ini¬ 
tial configurations are not identically distributed. In that setting we find, under 
mild conditions, that the particle will still create a number of reflecting structures, 
which will cause its motion to be periodic (see corollary]^. 

In contrast to this we show, that if there are correlations between the scatterer’s 
initial orientations, i.e. these orientations are not independent, then this can have a 
noticeable impact on the particle’s motion. Specifically, we show that it is possible 


^As is described in [^, this self-avoiding motion is very different from what is used in most 
computer simulations, since it is deterministically generated and not an ad hoc modified random- 
walk. 
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Figure 1. A subset of the honeycomb lattice H = is shown. 

The particle’s initial position r = (0,0) and velocity v = (1,0) 
are indicated by a triangle and arrow in this paper, respectively. 
The lattice is shown with a random initial configuration, in which 
closed circles indicate left scatterers and open circles indicate right 
scatterers, respectively. 


to have a random configuration with the same distribution as in our original model 
in this paper, but with local correlations between scatters. Interestingly, these local 
correlations have a global effect on the particle’s dynamics, causing it to shift from 
being periodic to having an unbounded motion (see [^). As one might expect, 
although the particle’s motion is based on a deterministic rule, its dynamics is also 
very much dependent on the specific statistical properties of the model’s initial 
configuration. 

The paper is organized as follows. In section we introduce the basic model, 
we will consider throughout the majority of the paper. In section we describe 
the type of structures the particle can create and annihilate, as it moves through 
the lattice. We show how these two processes can interact, in such a way, as to 
cause the particle to have a periodic motion. Based on these processes, we show 
in section 1^ that, for p G (0,1), the particle in the model we consider will have a 
periodic trajectory with probability 1. 

In sectionj^we describe the particle’s trajectory for p = 0,1, where the medium 
is initially homogenous and the particle has an unbounded and self-avoiding type 
of motion, as is considered in [^. We then investigate the transition between the 
particle’s self-limiting and self-avoiding motion, as p approaches 0 and 1. In section 
l^we study the case in which the model’s initial configuration of scatterers is non- 
identically distributed, where we find that, under mild conditions, the particle still 
has a self-limiting motion. In section we introduce other random configurations, 
in which the scatterers are not independently chosen, to show the effect of local 
(short range) correlations between initial orientations. We show this can lead to 
non-periodic dynamics. Some final comments, including a number of open question, 
are given in section 

We note that although the main results of this paper are proven rigorously, the 
paper is written so that it can be followed without the need for the reader to work 
through the various proofs, since the main ideas and techniques used to prove the 
results in this paper, are described in words, prior to each result. 
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2. The Flipping Rotator Model 


In this section we describe the particular flipping rotator model, we will consider 
throughout this paper. The lattice, over which the particle moves, is the honeycomb 
lattice H = (IHI,B), with sites El and bonds B. This lattice consists of regular 
hexagons, with sides of unit length, so that each lattice site has three nearest 
neighbors with which it shares a lattice bond of length 1 (see figure . 

Our main object of study in this paper is the motion of a single particle, as it 
moves from an initial position on the lattice, over the lattice, along the bonds of 
H, from lattice site to lattice site. By way of notation, we let r{t) G denote the 
position and v(t) € denote the velocity of the particle at time t > 0, where the 
particle is assumed to move with constant unit speed. Moreover, we let r = r(0) 
and V = v(0) denote the particle’s initial position and initial velocity, respectively. 

For simplicity, we also restrict ourselves to discrete time steps t = 0,1,2,..., 
so that the particle is at some lattice site at each time t. The particle’s trajectory 
is then the sequence of positions {r(t)}t>o C H. Since the velocity of the particle 
does not exist at the moment it is scattered, we let v(t) denote the velocity of the 
particle immediately after each time step t > 0. 

At each lattice site h G H, we assume that there is a scatterer, which rotates 
the velocity of the incoming particle, either to its left or to its right, by an angle 
of 9 = ±7r/3, depending on the scatterer’s orientation. This is shown in figure]^ 
where we use the convention, here and throughout the paper, that a closed circle 
denotes a left rotator and an open circle denotes a right rotator, respectively. The 
assumption, that there is a scatterer at each lattice site of H, ensures that the 
particle remains on the lattice for alH > 0. 

Note that each scatterer is, initially, either a left scatterer or a right scatterer, 
respectively, i.e. is oriented either to the left or to the right. With this in mind, we 
let C = C(0) denote the initial configuration of scatterers and let C{t) denote the 
configuration of scatterers on the lattice at time t > 0. The configuration C(t) for 
t > 0 consists of the collection of all the individual orientations of each scatter on 
the lattice at time t. For each lattice site h G El we let 


C{t, h) G {—1,1} for h G El, t> 0, 


denote the orientation of the scatterer at the site h at time t. The orientation 
C{t, h) = — 1 indicates, that at time t the scatterer at lattice site h is a left scatterer, 
whereas the orientation C{t, h) = 1 indicates, that the scatterer is a right scatterer. 
Furthermore, we let C'(h) = (7(0,11) denote the initial orientation of the scatterer 
at h G El at time t = 0. 

Suppose the particle has initial position r and initial velocity v. For an initial 
configuration (7, we call I = (r, v, (7) an initial condition. For an initial condition 
/, the particle’s deterministic equations of motion are given by 


r{t + l) = r{t)+v{t), 

v(f + 1) = i?[(7(t,r(t + l))]v(t), 

—(7(t, h) if h = r(t + 1) 

C(t,h) otherwise 


(7(t+l,h) = 


( 1 ) 

( 2 ) 

(3) 


for t > 0. Equation Q gives the dynamics of the particle, describing its piecewise 
linear motion between successive scatterings. The rotation operator R : { — 1,1} -G 
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Figure 2. Upon arriving at a left (right) rotator, indicated by 
a closed (open) circle, the particle’s velocity is rotated to its left 
(right) by an angle oi 0 = ±7r/3. The shaded circles indicate the 
particle and the arrow its position before and after being scattered. 


in equation ([^, is the rotation matrix given by 


R[z\ 


cos(f 0 ) sin(fz) 
-sin(f 0 ) cos(fz) 


where z S {—1,1}, 


( 4 ) 


which describes how the particle’s velocity is rotated, when it arrives at a scatterer. 
Equation ^ describes the flipping motion of the scatterers. 

Given an initial condition /, the particle’s motion over the lattice is uniquely 
determined for all < > 0 by equations This leads us to the following 

definition, which describes the general type of LEG we consider in this paper. 


Definition 2.1. Let {H, I) denote the LLG with the initial condition I = (r, v, C) 
and equations of motion given by the equations We will call this the flipping 

rotator model on the honeycomb lattice with the initial condition I. 


In general, one can consider two types of initial configurations. The first are 
fixed or deterministic initial configurations, which are specific configurations typi¬ 
cally given by some deterministic rule. For example, the initial configuration C in 
definition |2.1[ although arbitrary, is a deterministic configuration, since we assume 
that it has been specified. These deterministic type of configurations are considered 
throughout [^. 

The second type of initial configurations we consider are random initial configura¬ 
tions, which are generated according to some probabilistic rule. A major difference 
between a deterministic configuration of scatterers and one that is randomly gen¬ 
erated is that, a deterministic configuration is a single unique configuration, while 
a randomly generated configuration is a realization of a random process, of which 
there are typically many realizations. 

In this paper, our main focus is on random initial configurations. The particu¬ 
lar type of initial configurations we consider here, will be generated based on the 
following probabilistic rule. Let {H;p) denote the flipping rotator model, in which 
each scatterer is initially a right scatterer with probability p, for p G [0,1]. We 
furthermore assume that, the initial orientation of one scatterer does not influence 
the initial orientation of any other. Hence, the collection of initial orientations C 
in the {H;p) model forms a collection of independent identically distributed (i.i.d.) 
random variables, in which the probability 

P[C(h) = l]=p for all h G H. 


( 5 ) 
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Since this model has a random initial configuration we may assume, without loss 
of generality, that the particle’s initial position and initial velocity are 

r = (0,0) and v = (1,0), (6) 

respectively, as is shown in figure]^ That is, the particle in the {H\p) model is at 
the origin of the lattice at time zero and moving to the right. 

Under these assumptions we define the following model, which is similar to the 
hipping rotator model described in definition |2.1| The main difference is that its 
initial configuration is randomly generated. 

Definition 2.2. For p G [0,1], let {H,p) denote the LLG with the random initial 
condition I = (r, v, C) given by equations and equations of motion given 

by equations 0-§. We call this the flipping rotator model on the honeycomb 
lattice with the random initial condition I. 


Although the particle’s dynamics are deterministic in the {H;p) model, the 
model’s initial configuration is random for each p G (0,1). Therefore, the particle’s 


motion is what is referred to as a deterministic walk in a random environment 18 
or, alternatively, in a random medium. Here, we refer to this type of motion as a 
pseudo-random walk, since it has aspects of both random and deterministic motion. 
An alternative, but equivalent, way of describing the particle’s motion in the {H;p) 
model is to say that, the particle’s trajectory is a random walk up to the point at 
which it returns to a lattice site it has already visited. Once the particle returns to 
a site it has visited, its motion is determined by equations 0-0. 

Despite the randomness introduced by the initial configuration in the {H;p) 
model, we find the particle’s trajectory, over the lattice, is qualitatively the same 
for any realization of this initial configuration. Specifically, for each p G (0,1), the 
particle’s trajectory in the {H;p) model will always be periodic. That is, although 
the random medium, through which the particle moves, will initially cause the par¬ 
ticle to have a random-like motion, ultimately the particle’s deterministic behavior 
will overcome this randomness and its motion will be periodic. 

By periodic we mean that, there is a tp < oo, such that the particle’s position 
satisfies r(t) = r{t -\- tp) for each t > 0, where tp is the particle’s period. The fact 
that the particle’s motion in the {H-,p) model is always periodic for p G (0,1), is 
one of the main results of this paper, which can be stated as follows. 

Theorem 2.3. (Self-Limiting Trajectories) For p G (0,1), the particle’s tra¬ 
jectory in the {H',p) model will be periodic with probability 1. 

The proof of theorem |2.3| is postponed until section The reason for this is, 
that we first need an understanding of the kinds of “structures” the particle can 
create, as it moves through the lattice and how these structures ultimately limit the 
particle’s motion. Here, the term structure refers to a collection of scatterers, which 
collectively have a specific effect on the particle’s motion. The type of structure, 
which is most important to the particle’s dynamics in this model, is what we will 
refer to as a reflecting structure. 

A reflecting structure is, in effect, a collection of scatterers that causes the parti¬ 


cle to reverse its trajectory back to its initial position (see definition 3.1). A number 


of reflecting structures can be seen in figure where a particular realization of the 
particle’s trajectory in the {H-,p) model is shown for p = 0.5, 0.7,0.9, respectively. 
In each case there are two reflecting structures iki and '^2 which, acting together. 
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Figure 3. A realization of the particle’s trajectory in the {H;p) 
model for p = 0.5,0.7, and 0.9. In each case the particle’s tra¬ 
jectory is periodic, which is due to the formation of two different 
reflecting structures and ^' 2 , between which the particle’s entire 
trajectory oscillates. The periods of these trajectories are respec¬ 
tively tp = 11516, 7944, and 22504. 


cause the particle to have a periodic trajectory and thereby “limit” the particle’s 
motion to a finite subset of the lattice (see theorem 3.7). 

In the following section we will show that, as the particle moves through the 
lattice it, acting as an “architect” of sorts, can either create, transform, or annihilate 
a reflecting structure. Moreover, each one of these processes is an integral part of 
particle’s dynamics in the {H\p) model and each one will be used to prove theorem 

ESI 

Before continuing, it should be pointed out that in other flipping rotator mod¬ 
els, where the particle moves on some lattice other than the honeycomb lattice, 
the particle will have a different motion, than that described in theorem |2.3[ For 
instance, on the triangular lattice it has been shown that, for any initial configura¬ 
tion of scatterers on the lattice, the particle will always propagate in a strip away 
from its initial position 13 . On the square lattice, in which every lattice site is 


occupied by a scatterer, it has been shown that the particle will have an unbounded 
non-periodic trajectory irrespective of the initial configuration of scatterers [^. 

Before finishing this section, it should be noted that the {H;p) model does 
not have a random initial configuration for all p-values. In particular, for p = 
0,1 the initial configuration is the deterministic configuration consisting of all left 
and all right scatterers, respectively, considered in the previous paper These 
configurations will be important in this paper and are discussed later in section 

Recall that one can think of the lattice, together with its scatterers, as a medium 
through which the particle moves. In this context, p can be used as a measure of 
how homogenous or inhomogeneous this medium is. With this in mind, one of our 
goals in this paper is to understand the difference between the particle’s dynamics 
when p £ (0,1) and when p = 0,1. That is, we are investigating how the degree 
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of homogeneity (or lack thereof), in a random lattice medium, affects the particle’s 
dynamics. In particular, we are interested in the dynamical transition the particle 
undergoes, as the lattice medium becomes more homogenous, i.e. as p approaches 
0 and 1. 


3. Creation, Transformation, and Annihilation of Reflecting 

Structures 


Before studying the {H-,p) model, we first consider the {H,T) model in which I 
is an arbitrary but fixed initial condition. The reason we do this is that, to fully 
understand the dynamics in the {H\p) model, we need first, to describe the type 
of structures that the particle can form in the {H, I) model. 

There are, in fact, two types of structures we will consider in this section. 
These are reflecting and semi-reflecting structures, respectively. Reflecting struc¬ 
tures cause the particle to return to its initial position via the same sequence of 
positions it took to reach the reflector, only in reverse. In contrast, a semi-reflecting 
structure may reverse the particle’s trajectory only part of the way back to its ini¬ 
tial position, at which point the particle begins moving along a different sequence 
of lattice sites. 

We begin by describing how the particle is able to create these structures. We 
then show how the particle can, first, transform one reflector into another reflector 
and, second, how the particle can annihilate a reflector. In the following section, 
section]^ the interplay between reflector creation, transformation, and annihilation 
will be used to prove theorem 2.3 (see section]^. 

To describe the structure of a reflector and semi-reflector, we need the following 
definitions. For the LLG let T[tifl 2 ] = {r(t) '-ti <t < t 2 }, so that T[tifl 2 ] 

denotes the trajectory of the particle from time ti to time t 2 - The particle in {H, I) 
is said to travel on a loop L = T[ti,t 2 ] from time ti to time t 2 , if i’(ti) = r{t 2 ) and 
r{t) ^ r(ti) for all ti < t < t 2 - In this case we call the position r(fi) = r(t 2 ) the 
base of L. 

The reflecting structures, that will be important for understanding the particle’s 
dynamics, are defined as follows. 


Definition 3.1. Suppose there are times 0 < ti < t* < ^2 in the LLG {H, I), where 

(a) Li = T[ti,t^,] and L 2 = T[t^,fl 2 ] are loops; 

(b) both Li and L 2 intersect r{ti -j- 1) and r{t 2 — I) exactly once; and 

(c) T[O,ti-I]nT[ti,t 2 ] = 0. 

Then we call the particle’s sequence of positions = T[ti,t 2 ] a reflecting structure 
of based at the point r(ti) = r{t 2 ). We say that the particle encounters the 

reflecting structure T = T[ti,t 2 ] at time ti. 

From definition |3.H the reflecting structure di consists of two loops, which are 
both based at the same point r(ti) = r(t 2 ). This point, which is also the base of 
the reflecting structure, is visited exactly three times by the particle as it moves 
through the reflector (cf. figure |^. 

These loops can have any (finite) shape or size on the lattice, with the following 
restrictions. First, the loops cannot intersect any of the particle’s trajectory, prior 
to the reflecting structure. Second, each of these loops must visit the two lattice 
sites r(ti -|-1) and r{t 2 — 1) adjacent to the reflecting structure’s base exactly once. 
This topological structure of a reflector is shown in figure a). An example of a 
reflecting structure on the honeycomb lattice H is shown in figure |^b). 
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(a) reflector topology 



(b) example of a reflecting structure 


Figure 4. The general topology of a reflecting structure is shown 
in (a). An example of a reflecting structure 'll and the particle’s 
return to the origin is shown in (b). In both cases, the base of the 
reflecting structure, at r(ti) = r(t 2 ), is indicated by a diamond. 
The reflector’s two loops are shown as dashed-dotted and dashed 
lines, respectively. Solid lines indicate the particle’s trajectory be¬ 
fore and after the reflecting structure, where the particle moves 
first away from and then returns to its initial position r. 


When the particle encounters a reflecting structure, it will move through the 
structure and exit it at the same point it entered, but in the opposite direction. 
From there, the particle will retrace its trajectory, prior to the reflecting structure, 
back to its initial position. In this sense, a reflecting structure acts as a “mir¬ 
ror”, which effectively reflects the particle back along its original trajectory. This 
physical description of a reflector’s effect on the particle summarizes the following 
proposition. 

Proposition 3.2. (Reflecting Property): — T[ti,t 2 ] is a reflecting structure 

of {H,I), then r{t 2 + t) = r{ti — t) for all 0 < t < ti. In particular, the particle 
returns to its initial position at time t = t 2 + ti. 

Before proving proposition |3.2[ we need the following notation. If 17 C 14, let 
Cq (t) denote the configuration of scatterers on the subset 17 C 44 of the honeycomb 
lattice at time t > 0. Moreover, let denote the configuration of scatterers on 

the subset S7 = 44 — S7, which is the complement of 17. 

The proof of this proposition is based on the following observation. Suppose 
the particle is scattered to the right (left) by a scatterer at some lattice site r(4), 
at time 4 > 0. If immediately after being scattered, the particle were, for some 
reason, to reverse its velocity and head back towards this scatter, it would find, 
upon returning, that the scatterer’s orientation had flipped from right (left) to left 
(right). This change in orientation would then cause the particle to be scattered 
back to the lattice site r(4 — 1), it had visited at time 4—1. 

That is, any event that causes the particle to be “reflected” back towards a 
lattice site it has previously visited, will cause the particle to reverse its entire 
trajectory back to its initial position. In this model, a mechanism that can cause 
such a reflection is a reflecting structure. This observation is the main argument 
used in the following proof of proposition |3.2[ 

Proof. Suppose ik = T[4i, 42 ] is a reflecting structure of (44, 4), where 4 = (r, v, C). 
Then r(42) = r(4i), since the particle is at the base of ik at times 4i and 42. 
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(^2 “1“ ^ “I" 1) h) — 


Moreover, parts (a) and (b) of definition 3.1 imply that v(t 2 ) = —v(ti — 1) and 
part (c) implies then that the configuration C^{t 2 ) = Cq,{ti — 1). 

Continuing by induction, suppose for some t S [0,ti —1), that r{t 2 +t) = 
v(t 2 + i) = —v(ti — t — 1), and C^{t 2 + t) = C^{ti — t — 1). Then equation Q 
implies that 

r(t 2 + t + 1) = r(ti -t)- v(ti -t-l) = r{ti -t-1). 

From equation ^ it then follows that 

—C^(ti — < — 1, h) if h = r(ti — t — 1) 

C^(ti —t — 1, h) otherwise, 

for each h e Hence, C^(t 2 + t + 1) = C^{tx — t). Additionally, we have 
v(t2 + ^ + 1) = i?[C(t2 + t, r(t2 + ^ + 1))] v(t2 + i) = 

—R[C{ti — t — 1, r(ti — t— 1))] v(ti — t — 1) = —v(ti — t — 2), 

since R[C{ti +t— l,r(ti +t— 1))] = +t — 2,r(ti +t — 1))]. 

By induction, it then follows that r{t 2 + t) = r(<i —t), v(t 2 +i) = —v(ti — f — 1), 
and C^{t 2 +t) = C^{ti — t — 1), for all 0 < t < ti — 1. For t = ti this implies that 

r{t 2 + h) = r{t 2 + ti - 1) + v(t 2 + ti - 1) = r(l) - v(0) = r(0), 
completing the proof. □ 


Proposition |3.2| states that after the particle exits a reflecting structure, its mo¬ 
tion will be on the same sequence of lattice sites, as it took to arrive at the structure, 
except in the opposite direction. Having established this “reflection property,” we 
give an example of the creation or formation of a reflecting structure. 

Example 1. (Reflector Formation) Consider the trajectory of the particle shown 
in figure\^h) , consisting of the solid, dashed, and dashed-dotted lines. Note that, by 


i tion\3.1 
i \3.^ thi 


definition\3 .1 
sition 


Ik = T[5, 27] is a reflecting structure so that, as guaranteed by propo- 
the particle’s position r(27 1) = r(5 — t) for 0 < t < 5. In particular. 


the particle returns to its initial position at time t = 32 by reversing its trajectory 
back to its initial position, as can be seen in the figure. 

There are three reasons why a reflecting structure returns the particle to its 
initial position, based on the proof of proposition |3.2[ The first is that the particle 
leaves the reflector from the same point it entered, but in the opposite direction. 
Next, because of the flipping motion of the scatterers, the particle then begins to 
retrace each one of its previous steps, back to its initial position. The third, is that 
Ik = T[ti,t 2 ], as a collection of lattice sites, does not overlap with the particle’s 
initial trajectory T[0,ti] before time ti. Therefore, the particle’s motion through 
the reflector does not affect the configuration of the scatterers along the particle’s 
return path, back to its initial position. 

If a reflecting structure = T[ti,t 2 ] does happen to intersect the particle’s 
trajectory prior to time ti, then the structure is no longer a reflecting structure in 


the sense of definition 3.1 (see part (c)). This structure is, however, still similar to a 


reflecting structure, both in terms of its topology and its effect on the particle, which 
we will describe below. We refer to such structures as semi-reflecting structures, 
which are defined as follows. 
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Definition 3.3. Suppose in the LLG that the sequence of positions $ = 

T[ti,t2] is a reflecting structure, with the exception that r(T) S T[ti,t2] for some 
0 < T < ti. Then we call $ a semi-reflecting structure of {H,I). 

The topology of a general semi-reflecting structure is shown in figure [^a) and 
an example of a semi-reflecting structure is shown in [^b). Similar to a reflecting 
structure, a semi-reflecting structure $ = T[ti,t 2 ] consists of two loops, both of 
which are based at the same point r(ti) = r{t 2 )- Additionally, both loops also visit 
the sites r{ti 1 ) and r{t 2 — 1) adjacent to their base exactly once. The difference, 
however, between a reflecting structure and a semi-reflecting structure is that, as 
the particle moves through a semi-reflector it passes through at least one lattice 
site it has already visited prior to encountering this structure. Therefore, part (c) 
of definition 13.11 does not hold for semi-reflectors. 

The physical effect a semi-reflecting structure has on the particle’s motion, is 
the following. When the particle exits a semi-reflecting structure, it will begin to 
reverse its trajectory back to its initial position, as if it had encountered a reflecting 
structure. However, since the semi-reflector intersects at least one of the lattice sites 
r(r) that the particle has previously visited, the particle, by moving through the 
semi-reflector, will have altered this scatterer’s orientation. The idea is that, when 
the particle arrives back at the lattice site r(T), the difference in the scatterer’s 
orientation at this site, will cause the particle to deviate from its original trajectory. 
In this way, the semi-reflecting structure will reverse the particle’s trajectory back 
towards its initial position but, before the particle arrives at its initial position, this 
process of reversal will stop. This is summarized in the following proposition. 

Proposition 3.4. (Semi-Reflecting Property): Let $ = T[ti,t2] be a semi- 
reflecting structure of If t < ti is the largest time such that r(T) G <i>, then 

r{ti —t)= r{t2 + 1 ) for all 0 < t < ti — t. 


To prove this proposition we first note that, if the particle encounters a semi- 
reflecting structure <&, then this structure does in fact act like a true reflecting 
structure up to the point when the particle returns to the lattice site r(r). Since 
can be considered, up to this point in time, to be a reflecting structure, proposition 


3.2 can be used to prove proposition 3.4 This is the main technique we use in the 


proof of this result. 


Proof. Suppose $ = T[ti,t 2 ] is a reflecting structure of {H,I) where / = (r,v, C). 
If r < ti is the largest time at which r(r) S $, let {H, fl-i-i) be the flipping rotator 
model with the initial condition R+i = (r(r -|- 1), v(r -|-1), C(t -|- 1)). Then $ is a 
reflecting structure of {H, It+i) and proposition implies that r{t 2 +t) = r{ti—t) 
for 0 < t < ti — T — 1. 

Moreover, for t = 1 — r — 1 the proof of proposition |3.2| implies that r{t 2 + ti — 
T — 1) = r(T -I- 1) and v(t 2 -f ti — r — 1) = — v(r). From equationit then follows 
that 

r(f2 + <1 - r) = r(r + 1 ) - v(r) = r(r), 

which completes the proof. □ 


An example of a semi-reflector and its effect on the particle is given in the 
following. 

Example 2. (Semi-Reflector Formation) Consider the trajectory of the par¬ 
ticle shown in figure \^b), consisting of the solid, dashed, and dash-dotted lines. 
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semi-reflector topology 



(b) an example of a semi-reflecting structure 


Figure 5. The topology of a semi-reflecting structure is shown 
in (a). An example of a semi-reflecting structure $ is shown in 
(b) and its affect on the particle’s trajectory. In both (a) and 
(b), the base of the semi-reflecting structure at r(ti) = r{t 2 ) is 
indicated by a diamond. The semi-reflector’s two loops are shown 
as dashed-dotted and dashed lines respectively. Solid lines indicate 
the particle’s trajectory before and after the reflecting structure, 
which demonstrate, in both figures (a) and (b), how the particle 
stops reversing its trajectory at the lattice site r(r). 


respectively. Here, d) = T[7,31] is a semi-reflecting structure and t = 2 is the 
largest time at which r(r) € $ and r < 7. Proposition S.f then implies that 
r(31 -f t) = r(7 — t) for 0 < t < 7 — t . However, as can be seen in the figure, 
r(38 — T -I- 1) ^ r(r — 1), so that beyond the time t = 38 — r, the particle no longer 
reverses its trajectory. Therefore, there is no guarantee that the particle will ever 
return to its initial position only that it will return to r(T). 

Returning to our investigation of reflecting structures, once the particle encoun¬ 
ters one of these structure, it will make its way back to its initial position. However, 
the particle’s motion does not stop at that lattice site, but continues on. Therefore, 
the particle may continue to encounter other reflecting or semi-reflecting structures. 

In the remainder of this section we study how a number of reflecting and semi- 
reflecting structures can interact to cause either: (i) the particle to have a periodic 
trajectory, (ii) annihilate each another, or (iii) form a combination of (i) and (ii). 

We begin by considering the case in which the particle first encounters the re¬ 
flecting structure iPi = T[ti,t 2 ], then the reflecting structure '^/2 = T[t 3 ,t 4 ] (cf. 
figure]^. In this case, according to proposition 3.2 after the particle leaves il '2 it 
will reverse its trajectory and pass again through iki, but in the opposite direction. 
The reason the particle reverses it trajectory through iki is that, the particle’s ini¬ 
tial trajectory through iki, causes the scatters of to change orientation. In this 
sense, the particle transforms dii as it passes through it. 

A consequence of this transformation is that, when the particle returns to this 
transformed version of Ti, it will reverse its original path through this structure. 
This second, but reversed sequence of positions through Ti is, by definition, the 
semi-reflecting structure d'l = Tjtsjtg] for some t^fle > t^. Based on proposition 
|3.2[ this structure, which we call the transform of has the property that 

r(t6 + t) = r{ti — t) for 0 < t < ti. 


(7) 



















SELF-LIMITING TRAJECTORIES OF A PARTICLE 


14 




(a) periodic trajectory (b) reflector annihilation 

Figure 6. A periodic trajectory of the particle, consisting of the 
two reflecting structures ihi and iE' 2 , is shown in (a). The parti¬ 
cle, with initial velocity to the right, first encounters the reflecting 
structure ihi, which reflects it towards 'h 2 - The particle, subse¬ 
quently, becomes trapped between the two reflectors. Here, the 
particle’s direction is not indicated, since the particle traverses each 
lattice bond of ifi and d >2 in both directions. In (b), the annihila¬ 
tion of the reflecting structure by the semi-reflecting structure 
<i>i is shown. Here, the particle stops reversing its trajectory at the 
site r(r), indicated by the dashed line. 

That is, a transform 'I'l of a reflecting structure ihi is again a reflecting structure 
in that, after the particle exits the transform dti, the particle will immediately 
reverse its trajectory back to its initial position in exactly the same way it did 
when it exited the original untransformed reflecting structure (see example]^. 
This notion, that a transformed reflecting structure is still a reflecting structure, is 
summarizes the following proposition. The proof of this proposition follows directly 
from equation which is itself a direct consequence of proposition |3.5[ 

Proposition 3.5. (Transform Property) If ^ = T[t 3 ,t 4 ] is a transform of the 
reflecting structure = T[ti, ^ 2 ], then r{t 4 + t) = r{ti — t) for 0 < t < ti. 

Proposition |3.5| describes an important general consequence of the particle’s mo¬ 
tion, when it moves through a reflecting structure di. As the particle passes through 
Ik, it modifies (flips) each of the structure’s scatterers so that, if the particle were 
to return, its motion through the structure would be reversed. In this way, the 
particle’s motion through dt causes to change into its transform ik, which we 
denote by writing ik i]/. This type of reflector transformation is illustrated in 
the following example. 

Example 3. (Reflector Transformation and Periodicity) Consider the par¬ 
ticle’s motion in the LLG shown in figure^^a). Here, the particle first encounters 
the reflecting structure iki = r[l,23], which causes the transformation iki 1 —>■ iki. 
The particle is then reflected back through its initial position and towards the second 
reflecting structure ^^>2 = T[24,46]. Once the particle passes through '^ 2 , causing 
vk 2 I—> 'k 2 , it returns to its initial position through the transform This passage 
through iki, therefore, causes the transformed reflector "i! i to transform hack to the 
original reflecting structure i.e. iki 1 —>■ iki. The particle then continues on 
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through the transform ^' 2 , causing this reflector to transform back to 2 o-nd the 
particle to arrive at its initial position at time tp = 92. 

The particle has then visited each lattice site ^' 2 , and each point in be¬ 

tween, an even number of times at time tp, implying C{tp) = C. Additionally, both 
r{tp) = r and ^{tp) = v. Since {r{tp),v{tp),C{tp)) = (r(0), v(0), (7(0)), then at 
time tp the system is in the same state it was at time 0. Therefore, the particle will 
repeat this same motion for all t > 0. That is, the particle’s motion in this LLG is 
periodic with period tp = 92. 


Thus, as demonstrated in this example, if the particle encounters two reflecting 
structures ihi and 'I' 2 , where the particle first encounters iPi, and does not return 
to before encountering ^> 2 , its motion will be periodic (see theorem 3.7 and 
figure [^a)). This is so, because the particle will first move through ifii and then 
'I' 2 , causing iPi i-A and 'I '2 ^ 2 - After this, the particle will reverse its 

motion through these two structures, moving first through the transform ifii and 
then through iiE' 2 , causing 'I'l i-A ifii and ^2 'I' 2 - The result is one period of the 

particle’s motion. The particle will then continue to repeatedly pass through both 
reflecting structures, in this same order. In this sense, these structures never lose 
their ability to limit the particle’s motion, i.e. they remain reflecting structures for 
all time, and are the reason for the particle’s periodic motion. 

This seems to suggest that, whenever a particle encounters two different reflecting 
structures, its trajectory will be periodic. This, however, is not always the case. 
Another possibility is that, before the second reflector is encountered, the first 
reflector is annihilated. This possibility is considered in the following example. 


Example 4. (Reflector Annihilation) Consider the particle’s motion in the 
{H,I) model, as is shown in figure^b). The particle’s trajectory in this figure is 
similar to the particle’s trajectory in figure^a), since the particle first encounters 
the reflecting structure = T[l,23], which reflects the particle back towards the 
second structure $1 = T[24, 50]. Here, $1 is not a reflecting structure but rather 
a semi-reflecting structure, since $1 nT[0,49] ^ 0. According to proposition 
there must then be lattice site r(T), at which the particle may stop reversing its 
trajectory, after encountering $ 1 . In this example, this is indeed the case, as can 
be seen in figure^b), where the particle’s trajectory following $1 is indicated by a 
dashed line. 

Based on this example, we observe that if the particle encounters first a reflector 
then second a semi-reflector, this is not sufficient to cause the particle’s trajectory 
to become periodic. This, by itself, should not be surprising, since a semi-reflector 
may only cause the particle to reverse its trajectory up to a certain point in time (cf. 
example^. What is perhaps surprising, though, is that after the semi-reflector $1 
intersects ihi, can no longer function as a reflector. That is, the semi-reflector 
$1 has destroyed iPi’s ability to reverse the particle’s trajectory back to its initial 
position. 

To see this, suppose that after the particle leaves ifi U $1, it encounters a sec¬ 
ond reflecting structure '^ 2 , as is shown in figure(^ If the reflecting structure ilii 
were, in fact, still acting as a reflector, after being annihilated, then the particle 
will become trapped by the reflectors ihi and 'i> 2 , as in example^ However, the 
particle does not become trapped but, instead, its trajectory deviates at the lattice 
site r(r) from its previous trajectory, (cf. figure^. Therefore, when the particle’s 
trajectory intersects an existing reflecting structure, the result is that the reflector’s 


3.4 
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Figure 7. The particle first encounters the reflecting structure 
which is annihilated by the semi-reflecting structure <i>i. The 
particle then encounters the reflecting structures ^2 then ^' 3 , 
which cause the particle’s motion to become periodic. The par¬ 
ticle’s direction is not indicated, since the particle traverses each 
lattice bond of its periodic trajectory in both directions. 


key property, of being able to reflect the particle, is lost. In this case we say that, 
the reflector has been annihilated. 

On the other hand, if, after the particle leaves iPi U $1 for the second time, 
it encounters the reflecting structure ^' 3 , as shown in figure^ then the particle’s 
trajectory does indeed become trapped between the reflectors ^^2 and 'I' 3 . That is, if 
the particle has annihilated a reflector, the particle can still continue creating more 
reflecting structures, as it moves through the lattice and then eventually become 
trapped between two different reflectors. 

In the following, we formally define the notion of a reflector annihilation. 

Definition 3.6. Suppose = T[ti,t 2 ] is a reflecting structure of {H,I). Then the 
particle is said to annihilate 'I' at time r > t 2 , if the particle: 

(a) does not encounter a reflecting structure at any time t S (^ 2 , 1 "], where 

(b) T is the first time that r(T) G d/. 

If the reflector if is annihilated at time r > 0, then it is no longer a reflecting 
structure beyond this point in time. This notion allows us to prove the following 
theorem, which can be summarized as follows. If there are two reflecting structures 
at any point in time in the {H, I) model, then the particle’s motion in this LLG 
must be periodic. 

Theorem 3.7. Suppose that in the {H,I) model, there are two reflecting structures 
ifi = T[ti,t 2 ] and ^^>2 = ^[^ 3 ,^ 4 ] at some time t > t 4 > t^ > t 2 > ti. Then the 
particle’s trajectory is periodic with period tp = 2(^4 -f ^3 — t 2 — ti)- 

The reason this theorem holds is that, the particle will always become trapped 
between any two reflectors ifi and if 2 it has created, but not annihilated. Specif¬ 
ically, after the particle encounters ifi, it will be reflected back through its initial 
position and will eventually encounter if 2 . This second reflector will then reflect the 
particle back towards ifi, which will in turn reflect the particle back towards if 2 . 
These reflections will continue for all time, in effect trapping the particle between 
these two structures. This interplay of reflecting structures is the main tool used 
to prove theorem |3.7[ 
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Proof. Suppose = T[ti,t 2 ] and '^2 — T\t 3 ,t 4 \ are both reflecting structures of 
(-ff, I) at some time r where I = (r, v, C) and r > ^4 > ^ 2 - Since '$2 is a reflecting 
structure, proposition 3.2 implies that r(t 3 — t) = r{t 4 + t) for t G [Ojts]. Thus, 
the particle first encounters ihi at time ti and then 1112 at time t 2 , after which it 
encounters the transform il/i of at time ^4 + ta — ^ 2 - The particle then continues 
on, returning to its initial position r at time ^4 + ^ 3 . 

After the particle arrives at r, the claim is that its position satisfies the relation 


i ‘(^4 + f-a T i) — r(t 2 + + t) for 0 < t < ^3 — t 2 — ti, ( 8 ) 

so that the particle encounters the transform 1112 = T[t 4 + 2 t 3 —t 2 —ti, 2 t 4 +t 3 —t 2 —ti] 
of ^2 at time ^4 + 2t3 — t 2 — ti. If this is the case, then proposition |3.5| implies that 
r( 2 (f 4 + t 3 — t 2 — ti+1)) = r(t 3 — t) for 0 < t < ta so that this holds, in particular, 
for 0 < t < <3 — ^2 — ii- Hence, 


r(2(t4 + t3 — t2 — ti)) — r{t3 + ^ 2 ) — r, (9) 

v(2(t4 + t3-t2- h)) = v(f3 + t2) = V. (10) 


To verify note that proposition |3.2| implies that r{t 2 + ti) = r{t 4 + ta), 
v(t 2 +ii) = v(t 4 Tta), and C^^^{t 2 Ph) = C^^^{t 4 Tta). This last equation 
follows from the fact that the particle has visited each lattice site, except those of 
the reflecting structures an even number of times, by the times ^2 + and ^4 + fa- 
Given that r(t2 + G + t), r(t 4 + ta + t) ^ ^*1 U 1112 for 0 < t < ta — ^2 — G, then ^ 
holds, which in turn implies tha t ([^-( 10 ) hold. 

Note that each lattice site of ihi U 'I '2 has been visited an even number of times 
at time 2 (t 4 + ta — t 2 — G)- Moreover, the particle has passed through dii, ilii, 
'I' 2 , and ^2 exactly once by this time, so that each lattice site of U has 
been visited an even number of times at time 2(^4 + ta — ^2 ~ G)- Therefore, 
C(2{t4 + t 3 — t 2 — ti)) = C, which implies, together with equations that 

the particle has a periodic trajectory with period tp = 2(^4 + ta — ^2 — G)- HI 


This proof of theorem 3.7 explains why the particle in the (iJ, I) model becomes 
trapped between two reflecting structures. What is fundamental to this process 
is, that the scattering rule in this LLG is deterministic. That is, a right (left) 
rotator always rotates the particle’s velocity to its right (left) and then flips its 
orientation. If this were to happen, instead, with some probability p < 1, then any 
reflecting structure di could only have the reflecting property given in proposition 
|3.2| with some probability less than 1. In this case, the probability that the particle 
remains trapped between two reflectors will decrease to zero, as t —^ 00 , in which 
case theorem would not hold. Thus, the deterministic nature of the scattering 
rule is essential for the particle’s periodic dynamics in the {H ; p) model. 

Returning to the result of theorem |3.7[ once there are two reflecting structures 
in a given (77, 1) model, the particle’s dynamics is guaranteed to be periodic. Once 
periodic, the particle cannot visit any new lattice sites and, therefore, cannot en¬ 
counter any new reflecting structures. This limits the number of reflecting struc¬ 
tures, that can exists at any time in a flipping rotator model, to two. This argument 
proves the following corollary, which follows directly from theorem |3.7[ 


Corollary 1. At any time t > 0, in the LLG {H,L), there cannot be more then 2 
reflecting structures . 
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Although two is the maximal number of reflectors that can exist in the {H;p) 
model at any point in time, there is no limit to the number of different reflectors 
that a particle can encounter. For instance, in example the particle encounters 
three reflectors before becoming periodic. In fact, the reflector is possible that each 
time the particle creates a reflector, it is subsequently annihilated by the particle. 
In this way, the particle can encounter arbitrarily many reflecting structures as it 
moves through the lattice. 

In the following section, we will consider whether, in the {H;p) model, the 
particle is likely to encounter a reflecting structure and whether or not this structure 
is likely to be annihilated. What we find is that, ii p G (0,1), the particle cannot 
indefinitely avoid encountering a reflecting structure iki. This will be the first step 
in showing that the particle, for these values of p, must have a periodic motion. 
The second step is to show, that it is also unlikely that the reflecting structure 
is annihilated before the particle encounters a second reflecting structure il' 2 - 


4. Self-Limiting Motion in the { H ; p ) Model 


In this section, our goal is to give a proof of theorem |2.3[ which states that for 
p G (0,1), the particle’s trajectory in the {H;p) model is periodic. The proof of this 
result relies to a large degree on the results found in section specifically on the 
notions of reflectors, semi-reflectors, and transforms. Once a proof of this theorem 
has been given, we will study how the particle’s period depends on the parameter 
p. This will prompt us, in section to investigate how the particle’s dynamics 
abruptly changes as p approaches 0 and 1. 


To prove theorem 2.3 we will first show that the flipping rotator model {H,I) 


is invertible, i.e. its equations of motion are time-reversible. This time-reversibility 
will then be used to prove that the particle in the {H, I) model, for any initial 
condition J, has a periodic trajectory if and only if it stays in a finite subset of the 
(infinite) lattice H for all time. 

To show that the particle’s equations of motion (§-(§ can be reversed, we note 
that these equations describe the particle’s motion in forward time. That is, given 
r(t), v(t), and C{t) we can compute each of these quantities at time t -|- I. 

In the following proposition, these quantities in reverse time are shown to exist, 
i.e. given r(t), v(t), C{t) these quantities can be found at time t—1, (see equations 
@-([T3l)). The fact that these equations exist implies that, the particle’s motion 
is time-reversible. Therefore, it is possible to recover the past states of the [H, I) 
model based only on its present state. This summarizes the following proposition. 

Proposition 4.1. (Time-Reversal) For the initial condition I = (r, v,^), the 
particle’s time-reversed equations of motion in the {FI, I) model are given by 

r(t- 1) = r{t) - R[C{t,r{t))]v{t), (II) 

v{t-l) = R[C{t,r{t))]v{t), (12) 

-C{t, h) if r{t) 

(^(t, h) otherwise, 

for t > 1. 


C{t-l,h) = 


(13) 


The proof of this proposition is based on the following observation. In the {H, I) 
model, if one knows the particle’s velocity v(t) at time t, then there are only two 
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possibilities for what the particle’s velocity v(t — 1) at time t—1 could have been. If 
the particle encounters a right rotator at time at time t, then the particle’s velocity 
v(t — 1) will be one of these two possibilities. If the particle encounters a left 
rotator at time t then v{t — I) will be the other. Since it is possible to uniquely 
recover v{t — I), based on the type of scatterer the particle encounters at time t, 
it is possible to uniquely determine the particles position r{t — I) at time t — 1. 
Therefore, it is possible not only to know the particle’s future trajectory but also 
its past. A proof of proposition |4.I|is the following. 


Proof. Let I be the initial condition I = (r, v, C) and let r > 0. Because of the 
flipping motion of the scatterers, i.e. equation ([^, we have 

i?[C7(r + l,r(T + I))] = i?[- C'(T,r(r + 1))] = [^(t, r(T + I))]. 

Using this together with equation ([^, we have 

v(t) = R~^ [C(r, r(T + I))]v(r + I) = R[C{t + 1, r(T + 1))] v(t + 1). (14) 

Similarly, equation Q implies that 

r(T) = r(r + 1) — v(t). (15) 


Moreover, for each h G El we have 


C(r,h) 


-C(r + l,h) 

^(r + l,!!) 


if h = r(r + 1) 
otherwise. 


(16) 


since only the scatterer at h = r(r + 1) changes orientation at time r + 1. 

Replacing t by t — 1 in equations (14|-(16), yields equations Hence, 

equations (14)-(16l give the particle’s time-reversed equations of motion for t > 1 
completing the proof. □ 


Continuing towards our goal to prove theorem |2.3[ we will need to investigate 
what happens in the flipping rotator model if the particle’s trajectory is bounded. 
That is, suppose the particle in {H, I) remains in the subset oi H for all time. If 
n is bounded, i.e. finite, there must be two times 0 < U < t 2 at which r(ti) = r{t 2 ), 
v(ti) = v(t 2 ), and C{ti) = C'(< 2 )- The reason for this is that, because of the discrete 
nature of the lattice, the particle can only assume a finite number of positions and 
velocities on fl. Moreover, there are only a finite number of scattering configurations 
possible on fl. Therefore, at some time t 2 , the particle’s position, velocity, and the 
configuration of scatterers on the lattice must be the same as at some previous 
point in time ti < t 2 . To summarize, if the particle’s motion is bounded, its motion 
is either periodic or eventually periodic. 

Formally, the particle’s motion in (H, I) is said to be eventually periodic with 
period tp < oo, if there is a T > 0, such that r(t) = r{t + tp) for all t > T. 
Importantly, if T = 0, we do not consider the particle’s motion to be eventually 
periodic, since it is then periodic with period tp. 

If the particle, in the LLG {H, I), has either a periodic or eventually periodic mo¬ 
tion, its trajectory must be bounded. Conversely, using the time-reversal property, 
i.e. proposition |4.1[ it follows that, if the particle has a bounded trajectory, then 
its motion must be periodic (see for instance 20 , page 2). Hence, the particle’s 


motion in the LLG {H, I) cannot be eventually periodic. 

This argument is enough to prove the following proposition, which states that 
bounded motion and periodic motion are identical in the {H, I) model. That is, the 
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particle remains in a finite subset of the honeycomb lattice only if it has a periodic 
trajectory. 

Proposition 4.2. (Boundedness Property) The particle’s trajectory in {H,I) 
is periodic if and only if the particle’s motion is bounded. 

With proposition |4.2| in place, we are now in a position to prove theorem |2.3[ 
which says that the particle in the {H;p) model will, with probability 1, have a 
periodic motion if p € (0,1). Our strategy for proving this, is to show that, if 
the particle does not have a periodic, but rather an unbounded trajectory, then 
the particle will eventually encounter two reflecting structures, neither of which is 
annihilated. However, this will imply, via theorem |3.7[ that the particle’s trajectory 
is in fact periodic, contradicting our original assumption that it had an unbounded 
trajectory. This contradiction will be enough to prove that we will always observe 
periodic motion in the (i7;p) model if p G (0,1). This summarizes the proof we 
now give of theorem |2.3| 


Proof. Let Hi C iL be a single hexagon of the hexagonal lattice that contains the 
origin. For each n > 1, let be a collection of hexagons, such that C 

C iL, where there are exactly five hexagons between the boundary of 
H„ and the boundary i9H„_|_i of as is shown in figure]^ Let 'Ll and 'I '2 be 

the reflecting structures shown in this figure, where '^2 is any one of the reflectors 
labeled 4 > 2 . 

For the reflectors 'Ll = T\ti,t 2 ] and '^2 = T[t^,t 4 \, let )Li = T[ti — 5,^2 + 5] 
and _^i = T\t^ — 6,^4 + 6 ]. We call _^i and ^2 i-ii® extensions of the reflectors ^^i 
and 412 , which are those parts of the particle’s trajectory, containing ^^i and 4>2, 
respectively, in the set H„_|_i —H„. In particular, the lattice sites r(ti —5) = r(t 2 + 5) 
of )Li and r(t 3 — 6 ) = r(t 4 + 6 ) of )L 2 are on the boundary of i9H„. We call these 
the base of )Li and 4 ( 2 , respectively. (These are labeled r(Ti) and r(r 2 ) in figure]^ 
respectively.) For convenience, we say that any extended reflecting structure is a 
version of 4 [^ 2 ; it i® tti® same as '^2 to translation, reflection, and (or) reflection 
(see figure where a number of versions of _^2 shown). 

For p G (0,1), suppose the particle exits H„ for the first time at time ti > 0. The 
particle will then encounter the extended reflecting structure dfi = T\ti — 5, t 2 + 5], 
as is shown in figure]^ with probability P()Li) = p’'(l — pY > 0 , where r and £ are 
the number of left and right scatterers, respectively, that are initially in This 
is shown in figurewhere r(ri) = r(ti — 5) = r(t 2 + 5). 

According to proposition |3.2[ once the particle encounters 'Ll, it returns to its 
initial position at time ti +£ 2 - If the particle then remains in for all t>ti +t 2 , 
its motion is bounded and therefore periodic by proposition |T^ If the particle does 
not remain in, but exits H„ a second time, at time r 2 > ti, there are two cases. 
Either r(Ti) = r(T 2 ) or r(ri) ^ r(r 2 ). 

Case I: Suppose r(ri) = r(T 2 ). Then v(ri) = v(t 2 ) and moreover the configu¬ 
ration C'^Yun„(H) = C'^ 7 un„('r 2 )- Hence, r(Ti +t) = r(r 2 -f t) for t G [0,ti - ti], 
where r(r 2 -l-ti — Ti) = r(ti) is the base of the nonextended reflecting structure 'Ll. 
Therefore, at time T 2 -I- ti — Ti, the particle encounters the transfo rmed reflector 
'Ll = T[r 2 -l-ti — n, r 2 -l-t 2 — Ti] of 'Ll, which according to proposition 3.5 will cause 
the particle to reverse its trajectory back to its initial position. 

In particular, at time tp = 2 (i 4 + t^ — t 2 — H), where £4 = T 2 + t 2 — Ti and 
ts = T2 + ti — Ti, we have that r{tp) = r, v(£p) = v, and C{tp) = C, implying that 
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C r^n+i described in the proof of theo- 
2.3 is shown, along with the reflecting structures iPi = T[ti,t 2 ] 


Figure 8. The region 
rem 


and "^>2 = T[t 3 , ^4], as well as their extensions = T[ti — 5, <2+5] 
and ^2 = ^[^3 — 6, <4 + 6]. A number of versions of ^2 8-'^® shown. 
These extensions are based at r(ri) and r(T2), respectively. 


the particle’s trajectory is periodic with period tp. This follows, using the same 
argument used to prove theorem |3.7| 

Case II: Suppose r(ri) ^ r(T2)- By construction, some version of the extended 
reflector 5;^2 = ^[^3 — 6,<4 + 6], based at r(r2), does not intersect any lattice site 


of (cf. figure [8 
]P(^2) = P’'(l -P) 


. The particle encounters this version of T2 with probability 
where r and £ are, in this case, the number of left and right 
scatterers that are initially in ^21 respectively. 

Since the extended reflecting structures and 3'2 overlap spatially, the 

reflector ilii is not annihilated before the particle forms the second reflector il'2. 
Hence, the particle must have a periodic trajectory, via theorem 3.7 However, this 


happens only with some positive probability P(di2) > 0. If this does not happen, 
then there are two possibilities. Either the particle remains in for all time, 

so that the particle’s motion is periodic, or the particle exits at some future 
point in time. If the latter is true, we can repeat our previous argument to show 
that, with some positive probability, the particle’s motion will be periodic. This 
can be done as follows. 

Let En be the event that if the particle exits twice, at two different lattice 
sites, it first encounters the extended reflecting structure C and later 

some version of ^2 + ^n+i, where ^4 n ^2 = 0- Then 

nEn) > P^l) X P(S'2) > 0, (17) 

since (1) the assumption that p S (0,1) implies P(3li) > 0 for i = 1, 2 and (2) the as¬ 
sumption that the scatterer’s orientations are independent implies that the particle 
encountering and 3(2 are independent events. This assumption of independence 
in the (77; p) model, furthermore, implies that {En} is a collection of independent 
events. Moreover, (17) implies that = ^{En) = 00 . The second Borel-Cantelli 
lemma (see [^, p. 201, for instance), therefore implies that P(limsup„_,,oo En) = 1 
so that, with probability 1, the particle encounters two reflecting structures, neither 
of which is annihilated. 
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(b) 

Figure 9. A realization of the particle’s trajectory in the (iJ, 1/2) 
model is shown in (a), along with the reflectors ihi and 'I'2- The 
particle’s sequence of displacements {|i‘(t)|}t<tp for tp = 10728 is 
shown in (b), for one period of its motion. The times on the t- 
axis at which the particle encounters and transforms the reflectors 
^-1 = r[944,998], ^'2 = T[3602,3704], = T[6308,6362], and 

'I'2 = T[8966,9068] are indicated. The thick black lines indicate 
the times at which the particle’s motion is a pseudo-random walk. 
Thin black lines indicate the particle’s subsequent reflections back 
to its initial position. 


Therefore, either the particle becomes trapped in some for all time and has 
a periodic trajectory by proposition |4.2| or, with probability 1, at some point in 
time, there must be two reflecting structures in the {H;p) model. In the latter 
case, theorem [2.3| implies that the particle’s trajectory is periodic. Since these two 
cases are the only possibilities, this completes the proof. □ 


Theorem 2.3 guarantees that the particle’s motion in the {H\p) model will always 
be periodic, if p € (0,1). To illustrate how this periodicity arises, we consider a 
realization of the particle’s trajectory in the {H\p) model, shown in figure where 
p = 1/2. Here, the two reflecting structures iki and 4>2, that cause the particle to 
have a periodic trajectory, are indicated. In figure [^b) the particle’s sequence of 
displacements {|i'(t)|}t<tp are shown, where tp = 10728. Additionally, the times at 
which the particle encounters and transforms iki and 'I'2 are shown, as well as the 
times at which and ^2 are transformed back into iki and '^> 2 , respectively. 

Recall that, once the particle encounters a reflector, its future motion is then 
fixed, up to the point in time when it returns to its initial position (see proposition 
3.2). However, before it encounters a reflector, its motion can be considered to 
be a pseudo-random walk, in the sense that its motion is partly random, due to 
the random initial configuration of the scatterers, and partly deterministic, due to 
the deterministic rules of motion. Hence, in figure the particle’s motion is a 
pseudo-random walk until time t=944, when it encounters After returning to 
its initial position at time t = 1942, the particle’s motion is again a pseudo-random 
walk until it reaches 4>2 at time t = 3602. 
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p-axis 


Figure 10. The average period tav { p ) of the particle in the { H ] p ) 
model is shown for p S (0,1). 


Once the particles encounters this second reflector ^ 2 ^ its motion, for all future 
times, is completely determined by its previous motion through the lattice (see 
theorem 3.7). Hence, each reflector causes the particle’s motion to abruptly transi¬ 
tion from a pseudo-random walk to a deterministic one. In fact, after the particle 
has encountered two reflectors, which have not been annihilated, as is the case in 
figure l^a), the particle’s motion has become completely deterministic, so that the 
randomness of the initial configuration no longer affects the particle’s motion. 

Now, having established that the particle has a periodic motion in the {H;p) 
model so long as p e (0,1), we will investigate, how the particle’s period tp depends 
on the parameter p. Our motivation for doing so is to understand, how the particle’s 
dynamics changes in the {H;p) model as p —>■ 0,1. 

Using numerical simulations, the particle’s average period tav{p), averaged over 
realizations of the models initial configurations, is plotted for p € (0,1) in figure [I^ 
What is perhaps the most striking feature of this graph is the fact, that it appears 
that 

lim tavip) = oo. 

p-^0,1 

This suggests, that at p = 0,1 the particle’s motion is non-periodic and is therefore 
unbounded, by proposition |4.2[ 

This raises two questions. First, what is the nature of this unbounded motion 
and second, how does the particle’s motion make the transition from the periodic 


motion we observe for p S (0,1), to the unbounded motion we observe in figure 10 
for p = 0,1. Both of these questions are considered in the following section. 


5. A Discontinuous Dynamic Transition to Self-Avoiding Behavior 

In this section we investigate how the motion of the particle in the {H ; p) model 
changes as p approaches 0 and 1. That is, we will investigate, how the particle’s 
dynamics changes as the scattering configuration, i.e. the medium through which 
the particle moves, becomes more homogenous, until at p = 0,1 the initial config¬ 
uration consists of all left or all right rotators, respectively. What we will find is, 
that the particle makes a dynamic transition from the self-limiting periodic behav¬ 
ior, observed in section]^ to both a non-periodic as well as a self-avoiding mode of 
motion. Specifically, for p = 0,1, the particle’s trajectory in the {H\p) model is a 
self-avoiding walk between returns to the particle’s initial position to which it will 







SELF-LIMITING TRAJECTORIES OF A PARTICLE 


24 


return an infinite number of times (see theorem 5.1 1 . Moreover, based on numerical 
simulations, the particle’s trajectory is no longer periodic but unbounded (cf. figure 


lOl. 


To describe the transition between these two types of motion, we begin by not¬ 
ing that for p € ( 0 , 1 ), the model’s initial configuration is randomly generated, 
whereas for p = 0,1 the model’s initial configuration is deterministically generated. 
Specihcally, for p = 0,1 the {H;p) model’s initial configuration consists of either 
all left scatterers or all right scatterers, respectively. We note that, since both of 
these configurations are homogeneous, they are ordered in contrast to the random 
configurations considered in the previous section, for p G (0,1). As a consequence, 
the particle in the {H]p) model for p = 0,1 has a very different type of motion 
when compared with the values p = ( 0 , 1 ). 

In order to describe the dynamics in these models, we will use the notion of a 
self-avoiding cycle (cf. [^). Suppose that in the LLG there are two times 

ti < t 2 , such that each position r{t) is distinct for ti < t < t 2 and r(ti) = r{t 2 )- If 
this is the case, we say that the particle moves on the self-avoiding denoted 

by 7 = T[ti,t 2 ] from time ti to t 2 - Moreover, we call the position r(ti) = r(t 2 ) the 
base of the cycle. 

To describe the difference between the particle’s motion in the {H-,p) model for 
p G (0,1) and p = 0,1, we need the following theorem, found in (theorem 3.2). 
This theorem states that the particle’s trajectory in the {H,p) model for p = 0,1, 
can be described in terms of cycles. Specihcally, the particle’s entire trajectory can 
be decomposed into a sequence of self-avoiding cycles, each of which is based at the 
origin. 


Theorem 5.1. (Cyclic Decomposition Theorem For p = 0,1 there is an 
infinite sequenee of times {ri}i>o in the {H]p) model, with tq = 0 , such that 

(a) 7 i = T[Ti-i,Ti\ is a self-avoiding cycle based at the origin for i > 0 ; and 

(b) each 7 ^ is symmetric with respect to the line x = 1 / 2 . 


Therefore, between returns to the origin, the particle in the {H;0) and {H;l) 
models moves through the lattice without crossing its trajectory. In stark contrast 
to this, what we observe in the {H-,p) model, for all p G (0,1), is that the particle 
moves through the lattice until it encounters a reflecting structure. Then particle 
reverses its entire trajectory returning to each of the lattice sites it has previously 
visited, until it again arrives at its initial position. 

A natural question then is, why is there such a difference in the particle’s motion 
in the {H;p) model for p G (0,1) and p S 0,1. One immediate answer is that, for 
each p G (0, 1), there is always a positive probability that the particle will encounter 
a reflecting structure at some point in time t > 0 , which follows directly from the 
proof of theorem |2.3| However, if p = 0,1 then the probability for this is zero. 
Hence, it is not possible, in either the {H-,0) or (i?; 1) model, for the particle to be 
“reflected” back to its initial position by a reflecting structure. This observation is 
stated as the following corollary of theorem |2.3| and theorem |5.1| 

Corollary 2. No reflecting structures are possible in the {H;p) model for p = 0,1. 


The physical reason that there are no reflecting structures for p = 0,1 in this 
model is that, for these p-values, the particle’s trajectory is self-avoiding away from 

^We note that, sometimes these self-avoiding cycles are referred to as self-avoiding polygons [221. 
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its initial position. Since the particle must cross its own trajectory, away from its 
initial position, in order to create a reflector, there can then be no such structures. 
This is the main idea used in the following proof we give of corollary 


Proof. For p = 0,1, suppose that the particle in the {H\p) model encounters a 
reflecting structure ih = T\ti,t 2 ]. Part (a) of theorem 5.1 then implies that, the first 
loop Li — T[ti,tf\ of T (see definitionmust contain the origin, since r(ti) and 
r(t») are the same lattice site. The fact that r(0) is the origin and T[0, ti —l]n4' = 0, 
then leads to a contradiction of definition |3.1| part (c). Hence, the particle can never 
encounter a reflecting structure in the {H\p) model, if p = 0,1. □ 


Based on corollary no reflecting structures are possible in {H]p) for p = 0,1. 
For p S (0,1), what we observe numerically is that, the particle’s trajectory in 
the {H;p) model only becomes periodic after encountering a number of reflecting 
structures. This leads us to conjecture that the probability of ever encountering a 
reflecting structure is in fact equal to 1 for p S (0,1). This is in contrast to what we 
observe in these numerical simulations for p = 0,1, where the particle’s trajectory 
appears to be unbounded and therefore non-periodic (cf. figure 7 in [^). 

To physically understand the transition between the periodic self-limiting motion 
and the self-avoiding motion observed in the {H\p) model, we consider the case in 
which p G (0,1) is slightly less than 1. Then the concentration of scatterers that 
are initially oriented to the left is 1 — p, which is then slightly greater than 0. Since 
the particle in the {H;p) model and the {H; 1) model have the same trajectory only 
until the particle in the {H\p) model encounters a left scatterer then, over short 
time intervals, these two models have the same dynamics, with a high probability. 

Consequently, the particle’s transient, or short term dynamics in the {H\p) 
model, for p close to 1, consists, with a high probability, of a number of self¬ 
avoiding cycles, since the particle’s dynamics in the {H, ; 1) model consist entirely 
of self-avoiding cycles. As p approaches 1, the number of self-avoiding cycles in¬ 
creases as the particle’s trajectory in the {H\p) and the (i7; 1) model become more 
similar until, in the limit, the particle’s entire trajectory can be decomposed into 
a sequence of these cycles. It is in this manner then that the particle makes a 
transition from the periodic self-limiting motion, described in theorem |2.3[ to the 
self-avoiding motion described in theorem |5.1| Moreover, by symmetry, the same 
holds as p approaches 0. 

As a final remark, we note that, although the self-avoiding motion observed in 
the (77; p) model for p = 0,1, is very different from the periodic motion observed 
for p S (0,1), there is one striking similarity between the two. 

For any p G [0,1] the particle in the {H;p) model will always return to its initial 
position an infinite number of times. For p G (0,1), the reason is, that the particle 
has a periodic motion, so every point of its trajectory is visited an infinite number 
of times. For p = 0,1, the reason for this behavior is more subtle and is described 
in detail in section 6. This recurrence property is formally stated in the following 
corollary, which follows immediately from theorem |2.3| and theorem |5.1[ 

Corollary 3. (Recurrence Property) For p G [0,1] the particle in the {H;p) 
model returns to its initial position an infinite number of times with probability 1. 


Before continuing to the next section, we recall that in the {H;p) model the initial 
orientations of the scatterers are, as a set, a collection of i.i.d. random variables 
(see definition 2.2). In the following sections we investigate the particle’s dynamics 
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if these orientations are not identically distributed (see section]^, and alternately, 
if there are correlations between these orientations, i.e. these orientations are not 
independent (see section]^. 


6. Non-Identically Distributed Configurations 

In the {H ; p) model, considered in sections [2]^ the model’s initial configuration 
of scatterers was assumed to be both independent and identically distributed. In 
this section we consider a generalization of this model, in which the model’s initial 
orientations C are not identically distributed. What we find is that, as long as 
the initial orientations remain independent from each other then, under some mild 
conditions, the particle in this modified version of the {H;p) model, will still have 
a periodic trajectory. 

To make this precise, let F = {/i,..., /„} where fi : [0,1] —>■ [0,1] for each 1 < 
i < n. For this family of functions F, we let {H;F{p)) denote the flipping rotator 
model with initial position and velocity r = (0,0) and v = (1,0), respectively. This 
model’s initial configuration C is the random initial configuration where, for each 
h G H, we assume that 

P[C(h) = 1] = fi{p) for some fi G F. (18) 

Furthermore, we suppose that the initial orientations of the scatterer are indepen¬ 
dent of each other, so that C is a collection of independent random variables. 

Observe that if n = 1, F = {/i} is a single function. In that case, the collection 
of initial orientations C in the {F[;F{p)) model form a collection of i.i.d. random 
variables. In fact, if fi{p) = p, then {H;F{p)) is the {H-,p) model considered in 
sections In contrast, if n > 1, the collection of initial orientations in the 

{F[; F{p)) model need not be identically distributed, since there is then more than 
one function in the set F. 

However, for any n > 1, the particle’s motion in the {F[;F{p)) model will be 
periodic as long as the model’s initial configuration is still randomly generated, i.e. 
if each function fi{p) 7 ^ 0 , 1 , so that none of the scatterers are initially oriented to 
the right or left, with probability 1 . 

This result can be summarized as follows. If, in a flipping rotator model, the 
initial orientations of the scatterers are generated both randomly and independently, 
then it does not matter whether each orientation is probabilistically generated in 
the same way, the result will be qualitatively the same: the particle will have a 
periodic trajectory. This is formally stated in the following corollary. 

Corollary 4. For p G [0,1], the particle in the (iJ;F(p)) model has a periodic 
trajectory, with probability 1, if fiip) 0,1, for each fi G F. 

The proof of this result is based on the fact that, even if the initial orientations of 
scatterers in a given flipping rotator model are not identically distributed, reflecting 
structures still exist with a positive probability throughout the lattice. Eventually, 
the particle can and will become trapped between two of these structures. 


Proof. Let 41 ^ and 4'2 be the reflecting structures described in the proof of the- 
Under the assumption that fi{p) 7 ^ 0,1 for each fi G F, it follows that 


orem 


2.3 


P(4'), P('I'i), P(d> 2 ) > 0. From the proof of theorem 2.3 it follows that the particle 
either becomes trapped in the region for some m > 1 or encounters, with prob¬ 
ability 1, two reflecting structures, neither of which is annihilated. Proposition |4.2| 
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Figure 11. A realization of the particle’s trajectory in the 
{H;F{p)) model, described in example for p = 0.2,0.5, and 
0.7. In each case, the particle’s trajectory is periodic, which for 
each value of p, is due to the formation of two reflecting structures 
ihi and ^2- The periods are tp = 41127,3748, and 12450, respec¬ 
tively. 


and theorem [3.7| together then imply, that the particle’s trajectory in the {H-,F{p)) 
model is periodic with probability 1. □ 

We note that the only condition needed for corollary to hold, is that F is a 
finite collection of functions. Then, as long as each fi{p) ^ 0,1, the particle in the 
{F[\F{p)) model will, with probability 1, have a periodic trajectory. That is, the 
result given in theorem |2.3| can be extended to the case in which, the model’s initial 
configuration of scatterers is not identically distributed, but only independent of 
each other. 

Additionally, we note that, according to corollary]^ the scatterer’s initial orien¬ 
tation can be distributed in nearly any way since there are no restrictions on the 
functions in F. As an example of this result, we consider the following. 

Example 5. Let F = {/i,/ 2 }, where 

Mp) = \ cos(^p) and / 2 (p) = 1 - ^ cos(|p). 

As can he seen in figure [7| each lattice site h G H is of the end of exactly one 
horizontal lattice bond. We let IHI“ denote the collection of lattice sites h G H, that 
are on the left hand side of a horizontal bond and those lattice sites h G IH that 
are on the right hand side of a horizontal bond, respectively. If 

f/i(p) */hGH+ 

\/2(p)7/hGH- 

then, according to corollary^ the particle in the {H;F(p)) model will he periodic 
with probability 1, for each p G [0,1). As an illustration of this, a realization of the 
particle’s trajectory in the {H]F{p)) model is shown for p = 0.2, 0.5, and 0.7, in 
figurein\ (cf figure^. 
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Admissible Hexagonal Configurations 



Non-Admissible Hexagonal Configurations 


Figure 12. The seven admissible hexagonal configurations and 
the six non-admissible hexagonal configurations are shown. Open 
and filled circles represent right and left rotators, respectively. 


Importantly, flip) / 2 (p) for each p G (0,1], so that away from p = 0, the 
model’s initial orientation of scatterers is not identically distributed. Despite this, 
the particle’s motion in the (iJ;F(p)) model will still be periodic, with probability 
1, for p G (0,1) , which follows from corollary^ 

7. Admissible Configurations and Local Correlations 


In sections we considered the iH\p) model, in which the initial orientation 
of scatterers was both independent and identically distributed. In section we 
considered the (iL;F(p)) model, in which the initial orientation of scatterers was 
independent but not identically distributed. In this section we consider another 
variation of the iH\p) model, in which we introduce correlations between the initial 
orientations of the model’s scatterers, so that these are no longer independent of 
each other. 

For this model, which we call the (iL; A(p)) model, we show that although the 
initial orientation of the scatterers is still assumed to have the same distribution 
as in the {H;p) model, i.e. each scatterer is a right scatterer with probability p, 
these correlations lead to an absence of reflecting structures, and consequently to 
a qualitative change in the particle’s dynamics. Hence, the assumption that the 
orientation of each scatterer is independently chosen, is essential to the particle’s 
dynamics, in the iH;p) model. 

The initial configuration of scatterers we will consider in the {H\Aip)) model 
will be admissible configurations, which were first introduced in . An admissible 
configuration (7 is a particular type of configuration on the honeycomb lattice H, 
in which each hexagon H of the lattice has a particular orientation of left and right 
scatterers. To define an admissible configuration we note the following. 

Since there are six lattice sites on any single hexagon of H, there are 2® = 64 
possible configurations of left and right scatterers on each hexagon. In fact, up 
to rotation and reflection there are only thirteen different possible configurations, 
which are labeled (1)-(13) in figure 12 These can be separated into two distinct 


classes (l)-(7) and (8)-(13), which are referred to as admissible and nonadmissible 
hexagonal eonfigurations, respectively, as shown in figure |12[ 

This distinction between admissible and nonadmissible configurations on a single 
hexagon can be extended to a configuration C of the entire honeycomb lattice in 
the following way (cf. [^). 
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Definition 7.1. Let C be a configuration of scatterers on the hexagonal lattice. 
We say that the configuration C is admissible if the restriction of C to each hexagon 
of H has an admissible hexagonal configuration (cf. figure [I^ . Otherwise, we say 
C is a nonadmissible configuration. 

Our interest in admissible configurations is due to the following result, found 
in 1^. It states that, if the {H, I) model has an admissible initial configuration then, 
between returns to its initial position, the particle’s trajectory is a self-avoiding 
walk. 


Theorem 7.2. (Self-Avoiding Motion Suppose {H,I) has the initial con¬ 
dition I = (r, V, C) where C is admissible. 

(a) If the particle has the sequence of times {rijigx o,t which r(ri) = r then each 
7 i = T[Ti-i,Ti] is a self-avoiding cycle based at r for each i 

(b) If T is the last time at which r(r) = r, then {r(t)}(>T- is a sequence of distinct 
positions. 


Based on theorem 7.2 the particle’s motion in {H, I) is affected by whether 


the model’s initial configuration is admissible or not. In the {II;p) model, we will 
see that, with probability I, the model’s initial configuration is nonadmissible for 
p G (0,1) but is admissible if p = 0,1 (see proposition 7.3). Therefore, we do not 
expect to see any self-avoiding behavior in the (i7;p) model if p G (0,1), but are 
guaranteed this type of behavior if p = 0,1. 

From a physical point of view, a configuration must have a certain degree of 
order to be admissible. Since a randomly generated configuration, in which each 
orientation is chosen independently, will not be ordered^ such configurations will 
not be admissible, e.g. the initial configuration in {II;p) for p G (0,1). On the 
other hand, the much more ordered configurations of all left or all right scatterers, 
respectively, are admissible, e.g. the initial configurations in (i7;p) for p = 0,1. 
This dichotomy regarding which p-values lead to admissible initial configurations 
and which lead to non-admissible initial configurations in the (77; p) model, is stated 
as the following proposition. 


Proposition 7.3. For p G [0,1], let V[Cp = A] be the probability that the initial 
configuration in the {II;p) model is admissible. Then 

(a) f'[Cp = a1] = 0 for each p G (0,1); and 

(b) F[Cp = A] = 1 for p = 0, 1. 


To prove this proposition, we first calculate the probability that any given 
hexagon of the lattice will begin with an admissible hexagonal configuration. What 
we show is that, this probability is not equal to 1 if p G (0,1) and is 1 if p = 0,1. 
The idea is that, if any single hexagon has a probability less than 1 of having an 
admissible configuration, the probability is zero that every hexagon of El will have 
an admissible configuration. This summarizes the method used in the following 
proof of proposition |7.3[ 


Proof. For p G [0,1], let P[C'p(77)] be the probability that a hexagon 77 C 77 in 
{II;p) has an initial configuration (7), as shown in figure 12 up to rotation and 
reflection, for 1 < 7 < 7. Then one can compute that 


P[Ci(77)] = (1 -p)®, P[Cp"(77)] = 3p2(l -p)4, P[C'®(77)] = 6p\l - p)\ 
P[C4(77)] = (1 -p)®, P[C®(77)] = 3p2(l -p)4, P[C®(77)] = 6p2(l -p)^, 
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and P[C'J(H)]=/. 

The probability P[C'p('H) = A], that a hexagon T-L initially has an admissible con¬ 
figuration, is then the sum 
7 

P[Cp('H) = A] = Y , ^[CpW] = l-6p + 2 V - 54p3 + 72p4 _ 54^5 


By minimizing this function over all p G [0,1], it follows that V[Cp{'H) = A] < 1 


for any p G (0,1). 
Borel-Cantelli lemma 
each p G (0,1). Hence, part 


Since there are infinitely many disjoint hexagons on H 
page 201) then implies that 
(a) of proposition 7.3 holds. 


see 


21 


Ppp = H] = 1 for p = 0,1, then part (6) holds as well. 


the 

P[C'p = A] = 0 for 
Furthermore, since 
□ 


Proposition |7.3| states that those initial configurations, which are randomly gen¬ 
erated for p G (0,1), are nonadmissible, whereas the initial configurations, which 
are deterministically generated for p = 0,1, are admissible. Moreover, if the parti¬ 
cle moves on the lattice H, which has an admissible configuration, the particle will 
never encounter a reflecting structure (cf. corollary . Thus, a particle can only 
encounter a reflector, if the lattice configuration is initially nonadmissible. This is 
summarized in the following table. 


p-values 

(I) configuration type 

(II) dynamics 

(III) reflectors 

PG (0,1) 

probabilistic & nonadmissible 

periodic 

yes 

p = 0,1 

deterministic & admissible 

unbounded 

no 


Table 1. The table shows the differences in the {H\p) model’s 
(/) type of initial configuration, {II) dynamics, and {HI) when 
reflecting structures can or cannot be created by the particle in 
this model, for different p-values, respectively. 


Part (I) of table seems to suggest that, by generating an initial configuration 
using some probabilistic rule, the result will always be a nonadmissible configura¬ 
tion. In fact, if we introduce local correlations between the (77; p) model’s initial 
orientations, we can probabilistically generate admissible configurations. Although 
these configurations will be disordered and random, the fact that they are admissi¬ 
ble will mean, via theorem |7.2[ that the particle’s motion through the lattice will be 
a sequence of self-avoiding cycles. In order to describe the admissible, but randomly 
generated, configurations, which we will consider in this section, we first note the 
following. 

The honeycomb lattice H consists of an infinite number hexagons, whose sides 
and vertices make up the lattice’s bonds and sites, respectively. Let S denote the 
set of hexagons on the honeycomb lattice, which are shaded in figure |13[ As can 
be seen in this figure, each lattice site of H belongs to exactly one hexagon H G S. 
If we specify the orientation of each scatterer on the hexagon H for all H G S, this 
will generate a configuration of scatterers on the entire honeycomb lattice H. The 
following model can then be defined. 

Let {II\A{p)) denote the flipping rotator model, in which the particle has the 
initial position r = (0,0) and velocity v = (1,0). The model’s initial configuration 
C is the random initial configuration where, for each p G [0,1] and H G S, all 
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Figure 13. The shaded hexagons, shown above, comprise the set 
S C H oi hexagons, used to generate the random admissible con¬ 
figurations in the {H]A{p)) model. 


scatterers of "H are oriented to the right with probability p and to the left with 
probability 1 — p, respectively. Furthermore, we assume that the initial orientation 
of the scatterers on any single hexagon of S does not influence the initial orientation 
of the scatterers on any other hexagon of 5, i.e. they are independently generated. 

Under these assumptions, the initial orientation of scatterers in the {H;A{p)) 
model are identically distributed, since each scatterer has the probability p of being 
a right scatterer and probability 1 — p of being a left scatterer. These orientations, 
however, are not independently chosen but are correlated^ since each scatterer on 
any hexagon of S has the same initial orientation in this model. 

A specific consequence of the correlations found in the {H\A{p)) model is that, 
any realization of the initial configuration in this model will automatically be an 
admissible configuration (cf. figure [T^ . Theorem 7.2 then implies that the particle 
in the (i7; A(p)) model will have a self-avoiding motion away from its initial position. 
This self-avoiding behavior is described in the following corollary of theorem |7.2[ 
which also states that reflecting structures cannot be formed in the {H;A{p)), for 
any value of p. 


Corollary 5. For any p S [0,1], there will be, with probability 1, a sequence of 
times {rijigx in the {H-,A(p)) model, such that r(ri) = r and the following hold: 

(a) Each = T[ri_i,Ti] is a self-avoiding cycle, based at r for each i G I. 

(b) If T is the last time r(r) = r, then {r(t)}(>,- is a sequence of distinct positions. 

(c) The particle never encounters a reflecting structure at any time t>0. 


The proof of this result is based on the fact, that the initial configuration in 
the {H\A{p)) model will always be admissible. Hence, the particle in this model 
will have a self-avoiding motion (see theorem 7.21. Because of this, the particle will 
never encounter a reflecting structure since, in order to form a reflector, the particle 
has to cross its own path away from its initial position. This is impossible, if the 
particle’s motion is self-avoiding. This summarizes the following proof of corollary 

[SI 


Proof. Since any realization of the initial configuration in the {H]A{p)) model is 
admissible for p G [0,1], (a) and (b) follow directly from theorem |7.2[ while (c) 
follows using the same argument as is given in the proof of corollary [2] □ 

As an illustration of corollary]^ a realization of the particle’s trajectory in the 
(H; A(l/2)) model is shown in figure [T^a) for t < 20000. Since the particle does 
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Figure 14. In (a) a realization of the particle’s trajectory in the 
{H-,A{l/2)) model is shown for t < 20000. In (b) the difference 
between the particle’s mean square displacement A(t) and AsA(t) 
in the and {H;A{l/2)) models are shown, respectively, 

for t < 3000. 


not return to the origin by the time t = 20000 in this realization, the particle’s 
entire trajectory, up to this point in time, is a self-avoiding walk, as guaranteed by 
corollary 

The major difference between the {H\p) model and the {H\ A{p)) model is that, 
in the latter, the scatterers on each hexagon % £ S are not independently chosen 
to be either right or left scatterers, but are chosen to be either all left scatterers or 
all right scatterers. Despite this difference, the initial configurations of scatterers 
in both these models are identically distributed and, in fact, these distributions are 
the same, i.e. in both models P[C(h) = 1] = p for any h S H. Additionally, we have 
that {H;p) = {H;A{p)) for p = 0,1. Therefore, the only difference between these 
two models are the local correlations found in the initial scattering configuration of 
the {H;A{p)) model. 

What we observe in the (i4; A{p)) model, via numerical experiments, is that the 
absence of reflecting structures leads to a lack of periodicity in the particle’s motion. 
This can be seen in figure [l^b), where the particles’s mean square displacement 
is plotted for p = 1/2 in both the (i4; A(p)) and the (i4;p) models. The particle’s 
mean square displacement is given by: 

A(f) =(|r(t)p) for t > 0 in the (H-p) model; and 

A^(t) =(|r^(t)p) for t > 0 in the {H\A{p)) model, 

where r(t) denotes the particle’s position in the {H;p) model and rA{t) denotes 
the particle’s position in the {H-,A{p)) model. The average (•) is taken over all 
realizations of the random initial condition / = (r,v,C), in each respective model. 
The norm | • | is the standard Euclidean distance in the plane. 

In figure [T^b), we see that the particle’s mean square displacement A(t) in the 
{H;l/2) model asymptotically approaches the value V = 2000 as a consequence 
of the particle’s periodic motion. Physically, the value V is the average distance 
the particle is away from its initial position during one period of its motion, in the 
{H; 1/2) model. In contrast, the particle’s mean square displacement Ayi(t) in the 
{H; A{l/2)) model, can be numerically approximated by the function A^(t) « 
which suggests that the particle in this model has an unbounded trajectory. 
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Thus, by introducing local correlations into the initial configuration of the {H ; p) 
model, it is possible to drastically change the particle’s motion from being periodic 
to what numerically appears to be a non-periodic super-diffusive motion for p = 1/2 
(see figure [T^ b)). This may seem somewhat surprising, since the correlations found 
in the {H;A(p)) model’s initial configuration are local, in that only scatterers on 
the same hexagon can be correlated. Yet, these correlations are sufficient to change 
the global properties of the particle’s dynamics. 


8. Conclusion 


In this paper we have considered the deterministic motion of a particle moving 
on the honeycomb lattice H, whose sites are randomly occupied by left and right 
flipping rotators. When each individual scatterer is initially oriented to the right, 
independent of the other scatterers with probability p G (0,1), we show that the 
particle will have a periodic trajectory, with probability 1. In contrast, ifp = 0,1 the 
particle’s motion is a sequence of self-avoiding cycles between returns to its initial 
position, which numerical simulations suggest to be unbounded and, therefore, non¬ 
periodic. 

When p G (0,1), the particle can, as it moves through the lattice, create a 
number of structures. These structures, which we refer to as reflecting structures 
(see section]^, cause the particle to reverse its entire trajectory back through every 
one of the lattice sites it has previously visited. Importantly, two of these reflectors 
can trap the particle, limiting its motion to a finite subset of the lattice and cause 
the particle to have a periodic trajectory (see theorem 3.71. 

These reflecting structures and their role in causing periodic dynamics, have been 
previously observed in this and other LLG models [8 11 . In this paper though we 
describe, how the particle not only creates these structures, but how the particle 
can also transform and annihilate them. By studying the interplay of these three 
processes of creation, transformation, and annihilation, we are also lead to intro¬ 
duce the concepts of semi-reflecting structures as well as transforms of reflecting 
structures (see section]^. The processes of reflector creation, transformation, and 
annihilation, as well as the notions of semi-reflecting structures and reflector trans¬ 
forms, are each important for the particle’s dynamics and are used to show, that 
the particle’s motion will be periodic for p G (0,1), with probability 1 (see section 


In contrast, when p = 0,1, as was considered in the particle, between returns 
to its initial position, has a self-avoiding motion. Since this self-avoiding motion 
also appears to be unbounded, the difference between the particle’s periodic self- 
limiting motion, which occurs for p G (0,1), and its self-avoiding motion, which 
happens for p = 0,1, is quite striking. Hence, as p approaches either 0 or 1, what 
we observe is that the particle undergoes a discontinuous dynamical transition from 
a periodic, self-limiting motion to a non-periodic self-avoiding mode of motion. 

By combining the results in this paper, regarding the particle’s self-limiting mo¬ 
tion, with the theory previously developed in [^, which describes the particle’s 
self-avoiding behavior, we are able to qualitatively describe the particle’s transi¬ 
tion between these two types of motion. For p close to 0 or 1, we find that, while 
the particle’s transient, or short term, dynamics consists of a finite number of self¬ 
avoiding cycles, its long term dynamics is ultimately periodic. As p approaches 0 
and I, i.e. when the initial orientations of scatterers become more homogenous, we 
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observe that the particle’s period becomes increasingly longer, so that, in the limit 
as p approaches either 0 or 1, the particle’s motion is no longer bounded, but is 
entirely composed of an infinite sequence self-avoiding cycles. 

Furthermore, we have shown that for p = 0,1, it is not possible for the particle to 
form reflecting structures, suggesting that an absence of reflecting structures leads 
to non-periodic behavior. Stated another way, as the medium through which the 
particle moves, becomes increasingly homogenous, the particle becomes less able 
to build reflectors until at p = 0,1, where the lattice is completely homogeneous, 
these structures can no longer exist. 

To better understand why random initial configurations allow the particle to 
form reflecting structures, we also consider in this paper a generalization of our 
original {H;p) model. In that model, which we refer to as {H;F{p)), the initial 
configurations of the scatterers are not identically distributed, but are still inde¬ 
pendently chosen. We find that even without being identically distributed, the 
particle’s trajectory will still be periodic (see section]^. In contrast, if the ini¬ 
tial orientations of the scatterers are not chosen independently, this can drastically 
change the particle’s dynamics, even if these orientations are identically distributed 
(see section]^. Hence, the particle’s motion in the LLG we consider in this pa¬ 
per, although deterministic, depends on the particular statistical properties of the 
model’s initial configuration. 

The findings in this paper also lead to a number of open questions. One is, 
whether admissible configurations, including the random admissible configurations 
introduced in this paper, always lead to non-periodic and therefore unbounded mo¬ 
tion. Numerically, we observe this, but whether this can be proved mathematically, 
is unknown. This is related to a second open question, regarding the geometric 
nature of the self-avoiding cycles found in those flipping rotator models with an 
admissible initial configuration. In each case, these cycles appear to have a fractal¬ 
like structure, but it is unknown to what extent this is, or is not, the case. Indeed, 
an understanding of the geometry of these cycles and their dependence on the sys¬ 
tem’s initial configuration could be important in a number of applications, since 
self-avoiding walks are used to model the growth of crystals, polymers, etc. 22 ■ 25 


Additionally, there are questions regarding the reflecting structures, that are 
found in the models considered here. For instance, one of these questions is whether, 
reflecting structures are the only mechanism that can cause periodic motion. In¬ 
deed, if this could be shown to be the case, then we have an affirmative answer to the 
open question, whether or not admissible configurations always lead to unbounded 
trajectories, via corollary (see section]^. 

As we plan to show in later papers, the notions of reflector creation, transfor¬ 
mation, and annihilation, are also relevant for other two-dimensional LLG models, 
that have yet to be considered in the literature. 
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